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Verfahren imd Anordniang zur Detek^on unci zur Mess^^nq 
der Phase von periodischen Bxosxqnalen 

5 Die Erfindung betrifft ein Verfahren und eine Anordnung zur 
echtzeitfahigen zuverlassigen Detektion und Messung der Phase 
periodischer physiologischer GroBen Oder Biosignale. 

Im Stand der Technik sind Verfahren bekannt, die zur Bestim- 
10 mung der Phase ein iiber den zeitlichen Verlauf des Biosignals 
gleitendes oder sequentiell angeordnetes Analysef enster 
verwenden. Auf den im Analysef enster liegenden Signalab- 
schnitt werden auf der Fourier-Transformation basierende 
Methoden oder deskriptive Statistiken angewandt ._So werden 
15 beispielsweise in der Perimetrie periodisch auf leuchtende 
Lichtmarken definierter Intensitat und hinreichend hoher 
Frequenz (uber etwa 4Hz) benutzt, urn die Funktionsf ahigkeit 
des visuellen Systems zu uberprufen. Fur den Funktionstest 
wird das Elektroenzephalogramm (EEG) erfasst und die Reizant- 
20 wort auf den visuellen Stimulus beztiglich der Amplituden und 
der Phase analysiert. Die Phase der Reizantwort ist einer der 
entscheidenden diagnostischen Parameter in der Funktions- 
diagnostik. 

25 Bei den bisherigen Verfahren ist nachteilig^ dass die statis- 
tische Unsicherheit der Detektion bzw. die Ungenauigkeit der 
Messung sehr hoch ist. Die Unsicherheit und die Ungenauigkeit 
ergeben sich aus der Signaltheorie als Folge von und im 
Zusammenhang mit der Lange des Analysef ensters . Die Theorie 

30 besagt^ dass mit abnehmender LSnge des Analysef ensters die 
statistische Unsicherheit und damit die Ungenauigkeit zuneh- 
men, was in der praktischen Signalanalyse auch hinreichend 
bewiesen und bekannt ist. Fur ein statistisch besseres Ergeb- 
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nis musste zunachst die Fensterlange vergroBert werden. Aus 
der Physiologie jedoch ist bekannt, dass sich die Phase rela- 
tiv schnell andern kann und diese Anderungen auch diagnos- 
tisch relevant sind. Bei einem langen Analysefenster geht die 
5 wertvolle Information uber die Phasenanderung verloren und 
die statistische Unsicherheit des Messergebnisses nimmt 
infolge der Anderungen nicht zwangslaufig ab. 

Der Erfindung liegt die Aufgabe zugrunde, ein Verfahren und 
10 eine Anordnung anzugeben, mit denen es mGglich ist, einen 
kausalen Phasengang in periodischen Biosignalen mit einer 
gegenuber herkommlichen Verfahren besseren ZuverlSssigkeit 
und hoheren Geschwindigkeit bei gleichzeitig reduziertem 
Rechenaufwand zu detektieren und zu messen. 

15 

Erf indungsgemalS gelingt die Losung der Aufgabe dadurch, dass 
periodische Biosignale entsprechend ihrem physikalischen und 
physiologischen Ursprung erfasst werden, dass parallel zum 
analysierten biologischen System ein Zustandsbeobachter 
20 aufgestellt wird, dass die Ausgangsgr5fien des biologischen 
Systems und des Beobachters mit Hilfe eines Kalman-Filters 
ausgewertet werden und zur Bestimmung der Phase genutzt 
werden. 

25 Bei dem erf indungsgemSBen Verfahren wird die Phase eines 
periodischen Biosignals ermittelt und fQr f unktionsdiagnosti- 
sche Zwecke genutzt. So ist beispielsweise eine veriangerte 
Phase im Vergleich mit der gesunder Versuchsobjekte ein wich- 
tiger Hinweis auf funktionelle Probleme des untersuchten 

30 biologischen Systems. 

Bei der erf indungsgemaBen Anordnung ist parallel zum unter- 
suchten biologischen System, das mit einem Zustandsmodell 
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nachgebildet wird, ein Zustandsbeobachter angeordnet^ der 
entsprechend dem Systemmodell die Zustandsgrofie Phase basie- 
rend auf dem Kalman-Filter schatzt. 

5 Von Vorteil ist hierbei, dass die Schatzung der Phase 
kontinuierlich stattfinden kann und kein gleitendes oder 
sequentiell angewandtes Analysef enster notwendig ist. Damit 
wird die Analyse der zeitlichen Phasenanderungen erst 
moglich. Im Gegensatz zum relativ komplizierten theoretischen 
10 Hintergrund dieses Phasenschatzers ist die praktische Umset- 
zung einfach, Im Vergleich zu herkommlichen Verfahren ben5- 
tigt sie wesentlich weniger Rechenleistung, so dass eine 
Phasenschatzung in Echtzeit mOglich ist. 

15 Die Erfindung wird im folgenden anhand der theoretischen 
Herleitung und eines Ausf iihrungsbeispieles naher erlautert. 
In den zugehorigen Zeichnungen zeigen: 

Fig. 1 ein Blockdiagramm eines Beobachterkonzeptes; 



20 



Fig. 2 



ein Systemmodell einer erf indungsgemSBen Anordnung; 



Fig. 3 



ein prinzipielle Darstellung des Zustandsbeobach- 
ters zur Messung der Phase in periodischen Biosig- 
nalen; 



25 



Fig. 4 



einen Verlauf einer geschatzten Phase fUr eine 
Harmonische der Frequenz 8Hz und. Phase 2rad fUr die 
statischen Kalman-Faktoren 2 und 20; 



30 



Fig. 5 



einen Verlauf einer geschatzten Phase fur eine 
verrauschte Harmonische der Frequenz 8HZ und Phase 
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2racl mit einem SNR von OdB (unten) und dynamischem 
Kalman-Faktor (oben) ; 

Fig. 6 eine Phasenschat zung des Signals wie in Fig. 5 mit 
5 statischen Kalman-Faktor en; 

Fig. 7 Ergebnisse der Phasenschat zung (rechts) an realen 
Signalen (links); 

10 Fig. 8 eine additive Oberlagerung einer Harmonischen einer 
Frequenz von 8Hz mit dem Rauschen und gezielter 
Verstimmung der Analysef requenz (oben) und den 
Phasenverlauf mit Anstieg (unten) . 

15 

Ein biologische System, das ein periodisches Biosignal produ- 
ziert bzw. auf ein periodisches Eingangssignal antwortet, ist 
in Fig. 1 als Zustandsmodell „reales System'" dargestellt. Die 
nachfolgend genannten Zustandsgleichungen (1) und (2) 
20 beschreiben dieses System ( Fettgedruckte groBe Buchstaben 
stehen fur Matrizen, kleine fUr Vektoren) : 

x(t) = A • x(t) + B • u(t); x(0) = x(to) ( 1 ) 

y(t) = Cx(t). (2) 

25 

Fur die weiteren Betrachtungen wird ein additives Signalmo- 
dell angenommen, welches eine harmonische Schwingung und 
weifies normalverteiltes Rauschen summiert: 

3 0 y(t) = y . sin(a) t + cp(t)) + (t) ( 3 ) 

Das Ziel besteht darin, ein Systemmodell zu konstruieren. 
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dessen Variable x(t) die Phase <p(t) des zu untersuchenden 
Signals y(t) reprasentiert . Die Phase kann nicht direkt 
gemessen warden, da sie Argument einer trigonometrischen 
Funktion ist. Daher wird eine Hilf skonstruktion benotigt, 
5 Eine solche Konstruktion ist ein Zustandsbeobachter, der 
parallel zum untersuchten System angeordnet ist. Der Beobach- 
ter schatzt die Zustandsvariable mit Hilfe der Minimierung 
einer Fehlerf unktion, welche die Ausgange des realen Systems 
und des Beobachters vergleicht. Auf diese Weise kann nach 
10 abgeschlossener Fehlerminimierung die Zustandsvariable Phase 
direkt gemessen werden. 

Das Blockdiagramm des Beobachterkonzeptes ist in Fig. 1 
dargestellt. Da x(t) nicht direkt gemessen werden kann, wird 
15 Xo(t) im Beobachter geschatzt. Die innere Schleife im 
Beobachter minimiert den Fehler von Ym(t) bezuglich y(t) mit 
Hilfe der Korrekturmatrix K. Fiir den Beobachter ergeben sich 
dann die Zustandsgleichungen (4) und (5): 



20 Xm (t) = A . (t) + B - u(t) + K . [y(t) - y ^ (t)} , ( 4 ) 

yM(t) = CxM(t). (5) 

Aus (4) und (5) f olgt : 

25 XM(t) = (A-KC).x^(t) + B.u(t) + K.y(t). (6) 



Es wird davon ausgegangen, dass beide Systeme unterschiedli- 
che Anf angsbedingungen haben. Daraus ergibt sich der Beobach- 
tungsf ehler : 

30 

e(t) = x(t)-x^,(t). (7) 
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Der Beobachtungsf ehler verschwindet iterativ mit Hilfe der 
Korrekturmatrix K, so dass 

e(t) = 0 for t->oo. (8) 

Die Dynamik und die Stabilitat der Schatzung konnen mit der 
Dif ferentialgleichung des Beobachtungsf ehlers (9) beschrieben 
werden : 



e(t) = x(t)-x^(t). (9) 

Durch Umstellung und weitere Zwischenschritte erhalt man: 

15 e(t) = (A-K-C).e(t). (10) 

Entsprechend dem Signalmodell (3) ist damit zu rechnen^ dass 
das untersuchte Signal durch Rauschen gestSrt ist. Urn den 
Einfluss des Rauschens zu reduzieren, wird ein Kalman-Filter 
20 eingesetzt. Das Rauschen beriicksichtigend wird das System 
durch folgende Zustandsgleichungen beschrieben: 

Systemzus tand : x(t) = A • x(t) + B • u(t) + M • r,(t) (11) 
Systemausgang : y(t) = C- x(t) + rp(t) ( 12 ) 

2 5 Beobacht er : (0 = A • x^^ (t) + B • u(t) + K(t) • |y(t) - C - x^ (t)] , (13) 

wobei 

K(t) ist die Korrekturmatrix, welche zu erreichen hat, 
30 dass e(t) = x(t)-Xj^(t)-^0, 

e(t) ist der Beobachtungsf ehler, 
rs(t)ist das Systemrauschen, und 
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rp(t)ist das Prozessrauschen. 

Urn die Herleitung zu vereinfachen, wird angenommen, dass die 
5 Rauschanteile breitbandige Gaufische Nullmittelprozesse mit 
bekannten Kovarianzen sind: 

cov (t^tj) = E{rp(t0 r;(t2)}= Rp(t,) .8(t, -t^) 
covjt,, t^) = E{r3(t,).r;(t2)}= R3(t,)-6(tj - 1^) 

10 Die Rauschanteile sind voneinander unabhangig, also 

cov,^Jt,,t2) = 0. (15) 

Fur eine konsistente SchStzung von x(t) muss die Fehlerleis- 
15 tung mit Hilfe der Matrix K(t) minimiert werden: 

E^^ii) . e(t)} = E{ef (t) + e^(t) + • - • + (t)}= f (K(t))=Min. , (16) 

Unter Annahme der stochastischen Beziehungen beziiglich der 
20 Kovarianzen wird entsprechend dem Kalman-Filter eine geeig- 
nete Korrekturmatrix K(t) abgeleitet: 

K(t) = cov,(t)C^R;'(t). (17) 

25 Die Formel fur die Fehlerkovarianz coVe(t) kann vom Kalman- 
Filter abgeleitet werden: 

coVc(t) = Acov,(t)+coVg(t)A^-coVg(t)C^R;'(t)Ccov,(t)+M-R^ (18) 



BESTATIGUNGSKOPIE 



wo 2004/047610 T/EP2003/013355 



5 



- 8 - 



Schatzung der Phase: 

Das untersuchte Signal wird entsprechend (3) aus der Summe 
einer Harmonischen und Rauschen modelliert: 



y(t) = y„,(t) + rp(t) = y • sin(cot + <p(t))+ T^it) ,<o = 2nf. (19) 

Die Phase ergibt sich aus der Dif f erentialgleichung : 

10 9(t) = -a.(p(t) + r3(t) a>0. (20) 

Daraus ergibt sich das Systenimodel in Figur 2- Die Phase kann 
nicht direkt gemessen werden. Daher liegt parallel zum System 
ein Beobachter, in dem ein direkter Zugriff auf die 

15 geschatzte Phase <pM(t) moglich ist. Allerdings der nichtline- 
are Anteil yni(t) in (19) ist fur das Beobachterkonzept 
ungunstig. Eine geeignete Linearisierung yi(cp(t),t) wird 
benotigt, um eine lineare Beziehung zwischen der Zustandsva- 
riablen (puit) und dem Ausgang yM(t) entsprechend (5) herzu- 

20 stellen. Basierend auf der Taylor-Linearisierung kann der 
Beobachter folgendermaBen formuliert werden: 

<Pm (t) = -a • (Pm (t) + K(t) • (y(t) - y^ (t)) (21) 

25 yM(t) = y,(<PM(t),t). (22) 

Entsprechend (21) und (22) wird der Beobachter modellierf, 
wie in Fig. 2 dargestellt. Im Ergebnis der Linearisierung im 
Arbeit spunkt cps, wird (22) mit (5) verkniipft. Daraus ergibt 
30 sich der Faktor welcher in (17) verwendet wird, um den 

Korrekturfaktor K(t) zu bestimmen: 
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K(t) = cov, (t) . y • cos(cot + tp^ii)) - R^^ (t). (23) 
Als Arbeitspunkt wird die zu bestimmende Phase gewahlt 
5 q>B(t) = q>M(0> (24) 

und die resultierende Dif f erentialgleichung fur die Phase ist 
9m (t) = -a • <Pm (0 + cov, (t) • y - cos(cot + (Pw (0) • [y(t) - y • sin((ot + (p^ (t))] • R"" (t). (25) 

10 

Entsprechend (25) kann der Phasenschatzer modelliert werden, 
wie in Fig. 3 dargestellt. 

Zur Phasenschatzung in y(t) muss die Fehlerkovarianz berech- 
15 net werden* Aus (18) folgt: 

cove (t) = -2a • cov, (t) - y ' • cos^ (<ot + cp^ (t)) • cov^ (t) • R"* (t) + R, (t). (26) 

Gleichung (26) ergibt eine einfache LOsung, falls hoherfre- 
20 quente Anteile in der Fehlerkovarianz nicht berucksichtigt 
werden. Basierend auf (27) 

cos^ (co t + (t)) = 1 / 2 + cos(2(ot + 2<Pm (O) (27) 
25 kann (26) vereinfacht werden zu: 

covc(t) = -2a • cov,(t) - 1 /2 • y' . Rp>(t) - cov^ (t) + R^Ct) . (28) 

Bei geeigneter Wahl des Parameters a in (28) werden hochfre- 
30 quente Anteile in Folge temporaler Integration unterdruckt, 
also es liegt ein Tief passverhalten vor, Unter Berucksichti- 
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gung des Tiefpasses kann (25) vereinfacht werden: 

ci>M (0 = -a - (t) + cov, (t) . y . cos(cDt + ip^ (t)) - y (t) • R^* (t), (29) 

5 Damit wird der Beobachter, dargestellt in Fig. 3, verein- 
facht. Das in (29) vorgeschlagene System kann insbesondere 
zur Phasenschatzung von Harmonischen im Rauschen eingesetzt 
werden. 

10 In Fig. 4 ist der Verlauf der geschatzten Phase fur eine 
Harmonische der Frequenz 8Hz und Phase 2rad fur unterschied- 
liche Kalman-Faktoren dargestellt. Die Kalman-Faktoren sind 
statisch und betragen 2 bzw. 20. Wie der Grafik entnommen 
werden kann, mit einem niedrigeren Kalman-Faktor wird die 

15 Schatzung langsamer. Statische Kalman-Faktoren miissen dort 
eingesetzt werden, wo der Zeitpunkt der Phasenanderung nicht 
bekannt ist. 

In Fig. 5 ist der Verlauf (unterer Teil) der geschatzten 
20 Phase fur eine verrauschte Harmonische der Frequenz 8Hz und 
Phase 2rad mit einem SNR (Signal-Rausch-Verhaltnis) von OdB 
und dynamischen Kalman-Faktor (oberer Teil) dargestellt. Ist 
der Zeitpunkt der Phasenanderung bekannt, so kann der Kalman- 
Faktor so konstruiert werden, dass zunSchst der Anderung 
25 gefolgt und anschlielSend die Varianz der Schatzung minimiert 
wird. 

Zum Vergleich ist in Fig. 6 die Phasenschatzung des ' selben 
Signals wie in Fig. 5 mit statischen Kalman-Faktoren darge- 
30 stellt. 
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In Fig. 7 sind die Ergebnisse einer Phasenschatzung (rechte 
Spalte der Grafik) an realen Signalen dargestellt. Zur Stimu- 
lation eines beispielhaft untersuchten visuellen Systems 
wurden Lichtpulsf olgen mit der Wiederholrate von 8 Pulsen pro 
5 Sekunde verwendet, wobei abwechselnd einer Ruhepause (Zeit 
von 0 bis 2 der Zeitverlaufe der Signale oben links und unten 
links) die Lichtstimulation (Zeit 2 bis 5) folgte. Darge- 
stellt ist der Verlauf des EEG (Elektroenzephalogramm, linke 
Spalte der Grafik) von zwei okzipitalen Positionen nach 16- 
10 facher reizbezogener Mittelung. Beide Zeitverlaufe zeigen 
einen deutlichen Sprung der Phase nach Beginn der Lichtstimu- 
lation. 

Die Phasenschatzung wird problematisch bei stark verrauschten 
Signalen. Generell gilt zwar, dass die Phase robuster gegen 
Storungen ist als die Amplituden, wie dies schlieBlich auch 
in der Inf ormationstechnik bekannt ist. Allerdings ist in 
diesem Grenzbereich zunachst die Frage nach dem Vorhandensein 
- also der Detektion - einer kausalen Phase zu klSren, erst 
dann ware die Phase zu schatzen. 

In Fig. 8 ist eine Harmonische der Frequenz 8Hz additiv dem 
Rauschen uberlagert beginnend bei t=4s, wobei das SNR -lOdB 
betragt (oberer Verlauf in der Grafik) . Im Zeitbereich ist 
25 die Harmonische beztlglich ihrer Amplituden nicht nachweisbar. 
Verwendet man den Phasenschatzer mit einer gezielten Verstim- 
mung, hier mit einer Frequenz von 7.8Hz^ also urn 0.2Hz 
weniger als die Frequenz der Harmonischen, so ergibt sich im 
Falle einer kausalen Phase ein Anstieg von 1.2rad/s (unterer 
30 Verlauf der Grafik) . Dieser Anstieg kann in Kombination mit 
einem Diskriminator direkt zur Detektion des Signals genutzt 
werden. 



15 



20 
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Bezugazeichenliste 

a - Systemparameter^ im Beobachtermodell wahlbar 

B, C, R - Matrizen im Zustandsmodell eines Systems 

cp(t) - Phase im Systemmodell, zu schatzende Grofie 

(PM(t) - Phase im Beobachtermodell, messbare Grolie 

rp(t) - Prozessrauschen 

rs(t) - Systemrauschen 

u(t) - Eingangsvariable eines Systems im Zustandsmodell 

x(t) - Zustandsvariable eines Systems im Zustandsmodell 

^(t) - Zustandsvariable des Beobachters im Zustandsmodell 

y(t) - Ausgangsvariable eines Systems im Zustandsmodell 

yM(t) - Ausgangsvariable des Beobachters im Zustandsmodell 
yi(<p(t)/t) - Linearisierungsoperator fur Phase 
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Paten tanoprflche 



10 
15 
20 
25 

2. 

30 



Verfahren zur Detektion und zur Messung der Phase von 
Antwortsignalen (y(t)) eines Biosystems^ folgende 
Schritte umfassend: 

a) Multiplikation des Antwortsignals (y(t)), dessen 
Phase (<p(t)) bestimmt werden soli, mit einem ersten 
Faktor; 

b) Multiplikation des aus Schritt a) erhaltenen Produkts 
mit einem zweiten Faktor, welcher durch eine trigono- 
metrische Funktion reprasentiert ist, deren Argument 
sich aus dem Produkt der Frequenz des untersuchten 
Antwortsignals mit der Zeit, addiert mit der gemesse- 
nen Phase, ergibt, wobei die Frequenz der trigonomet- 
rischen Analysef unktion der Frequenz entspricht, bei 
der die Phase ermittelt werden soil, oder um einen 
bekannten Betrag von dieser Frequenz abweicht, 

c) Multiplikation der gemessene Phase mit einem dritten 
Faktor (a) ; 

d) Dif f erenzbildung aus dem im Schritt b) erhaltenen 
Produkt und dem im Schritt c) erhaltenen Produkt; 

e) Integration der im Schritt d) erhaltenen Differenz 
Uber die Zeit, wobei das Ergebnis dieser Integration 
die zu ermittelnde Phase des Signals reprasentiert 

f) Wiederholung der Schritte a) bis e) , bis ein Abbruch- 
kriterium erfiillt ist. 

Verfahren nach Anspruch 1, dadurch gekennzeichnet , dass 
der erste Faktor zeitlich konstant oder veranderlich 
gewShlt wird. 
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3. Verfahren nach Anspruch 1 oder 2, dadurch qekennzeichnet ^ 
dass bei der Dif f erenzbildung im Schritt d) das Produkt 
aus Schritt c) von dem Produkt aus Schritt b) abgezogen 
wird. 

5 

4. Verfahren nach einem der Anspriiche 1 bis 3, dadurch 
gekennzeichnet^ dass das Antwortsignal (y{t) ) einem 
Zustandsbeobachter zugefiihrt wird, welcher die Verfah- 
rensschritte a) bis f) ausfuhrt, urn eine gesch^tzte Phase 

10 (q)M(t)) zu bestimmen, wobei das Verfahren abgebrochen 

wird, wenn das Beobachterausgangsignal (yM(t)) weniger 
als einen durch eine Fehlerf unktion (cove(t)) vorgegebe- 
nen Fehlerwert (e(t)) vom Antwortsignal (y{t)) abweicht 
und wobei nach Abbruch des Verfahrens die geschStzte 

15 Phase (cpnCt)) gleich der Phase (<p(t)) des Antwortsignals 

gesetzt wird. 

5. Verfahren nach Anspruch 4, dadurch gekennzeichnet , dass 
die Schritte a) bis d) im Beobachter nach der folgenden 

20 Formel ausgefuhrt werden: 

(0 = ' (0 + cov, (0 • y . cos(fij^ + ip^ (0) • y(0 • r;' (t) . 

6. Anordnung zur Detektion und zur Messung der Phase von 
Antwortsignalen eines Biosystems, dadurch gekennzeichnet , 

25 dass die Anordnung einen Zustandsbeobachter umfasst, dem 

parallel zum Biosystem das Antwortsignal des untersuchten 
Systems eingespeist wird und welcher die Verfahrens- 
schritte nach einem der Anspruche 1 bis 5 ausfiihrt. 

30 7. Anordnung nach Anspruch 6, dadurch gekennzeichnet, dass 
der Zustandsbeobachter einen Kalman-Filter umfasst, mit 
welchem Storsignale aus dem Antwortsignal gefiltert 
werden. 
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